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ABSTRACT 

The sound radiated by isotropic turbulence is computed using inertial range scaling 
expressions for the relevant two time and two point correlations. The result depends on 
whether the decay of Eulerian time correlations is dominated by large scale sweeping or 
by local straining: the straining hypothesis leads to an expression for total acoustic power 
given originally by Proudman, whereas the sweeping hypothesis leads to a more recent 
result due to Lilley. 
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I. Introduction 


The sound radiated by isotropic turbulence was first computed by Proudman 1 very 
shortly after the appearance of Lighthill’s quadrupole theory 2 of sound radiation from an 
arbitrary stochastic velocity field. For the total acoustic power p in the far field, Proudman 
found 

p = ap(-) 5 £ ( 1 ) 

c 

where c is the speed of sound and p the density in the far field, e is the dissipation rate of 
the isotropic turbulence, u is defined by 

u 2 = --K (2) 

where K is the turbulence kinetic energy, and the presumably universal proportionality 
constant a, the Proudman constant, was expressed in terms of the spatial correlation 
function of the turbulence. 

Proudman’s calculation was re-examined by Lilley 3 who noted that Proudman had in- 
troduced approximations which suppress dependence on time correlations. This important 
observation led to a modification of Proudman’s formula 


P = 


. / flT \ 4 / \ 5 

a p\—) (-) £ 


u 


c 


(3) 


where temporal properties appear through the frequency 0 which is defined so that the 
time correlation function depends only on the product Or, where r is the time difference. 
In Eq. (3), L is the integral scale defined by L = u 3 /e. By invoking a correction 4 to the 
quasinormal hypothesis suggested by Lighthill, Lilley also found that the proportionality 
constant a' depends on the flatness of the single point velocity probability distribution 
function. However, experimental evidence 5 indicates that this correction is not large; 
therefore, this effect will be neglected here. 

In evaluating the integrals required by Lighthill’s theory, Lilley substituted empir- 
ical formulas for the relevant correlation functions, choosing in particular the two time 
dependence to establish the far field sound spectrum found in numerical simulations by 
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Sarkar and Hussaini. 6 Lilley has argued 7 that this procedure is reasonable since the tem- 
poral covariance depends strongly on the large eddy contribution and therefore should be 
accessible to direct numerical simulation. 

In this letter, Lilley’s calculation is repeated using Kolmogorov inertial range forms 
appropriate to high Reynolds number turbulence. This procedure leads to expressions 
for acoustic power which contain only universal inertial range quantities. An important 
feature of Lilley’s analysis is the dependence of sound emission on both Eulerian two time 
and two point correlations. Whereas there is little doubt that the spatial correlations 
must follow the Kolmogorov laws < rt(x + r)rt(x) >~ r 2 / 3 or equivalently E(k) ~ & -5 / 3 , 
there are two plausible alternatives for the time correlations 8 : temporal decorrelation 
of the Eulerian velocity field can be caused either by sweeping by the largest energetic 
scales or by local straining. A definitive discussion of these decorrelation mechanisms 
and their physical basis is given by Kraichnan 9 ; briefly, under the straining hypothesis, 
temporal decorrelation is a local property of inertial range scales, but under the sweeping 
hypothesis, it is a nonlocal property of the energy containing range of scales. Although 
the most recent discussions 10,11,12 favor the sweeping hypothesis, we examine both cases 
and find that the local straining hypothesis leads to Proudman’s formula Eq. (1), but the 
sweeping hypothesis leads, after an appropriate substitution for 0, to Lilley’s formula Eq. 
(3). 

II. Evaluation of the sound radiated by isotropic turbulence 

For time stationary sources, Lilley’s approximate form 3 of Lighthill’s theory predicts 

P = / dr U ( r ) ( 4 ) 

where the integrand U is given by 

<9 4 

U(r) = -Qrj-i C < Ml ( X + + ' r ) 2 “l( X ^) 2 > - < Uj > 2 ] | T=0 

r \ 4 

= ^I W ( r ’ T ) |r=0 (5) 
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where u\ is the velocity component in the direction of the line connecting the source and 
observation point. The quantity hi of Eq. (5) is evaluated in wavenumber space by using 
the quasi-normal closure in the usual way to obtain 

W(k ,r) = ^/ dpdqQ(p,T)Q(q,T) (6) 

5 ik=p+q 

where Q( k, r) is the correlation function defined by 


< Ui(k,t + T)uj(k',t) > = Q(k,T)Pij(k)6(k + k') 


and 


Pij{ k) = 8ij - k l k l k 


Following Lilley, we set k — > 0 in Eq. (6), corresponding to the far field approximation 
Then 

(2tt) 3 p a 4 4 f 2 2 


P = 


87 r 2 c 5 dr 4 5 
Now write the two time correlation as 


- / 4vrp dp Q(p } t'Y 


( 7 ) 


Q(p,t) = Q{p)r{p,r ) 


(8) 


where 


A7rp 2 Q(p) = E(p) 


is the energy spectrum and r has the similarity form 


(9) 


r (p , r) = r{rd{jpy) 


( 10 ) 


with 


_ / Vp 


Whereas the decorrelation time is a local inertial range property under the straining 
hypothesis, the sweeping hypothesis makes the decorrelation depend on the sweeping ve- 
locity V, which is not an inertial range property, but an entirely independent property of 
the energy containing range. 


( 11 ) 


3 



Substituting Eqs. (8)-( 1 1 ) in Eq. (7) leads to 


P = 


1 p d 4 


5 c 5 dr 4 J 0 

Differentiating Eq. (12) under the integral sign, 

1 f°° 

P = P 4 / d P p~ 2 E(p) 2 K 

5 c 5 J 0 


dp p 2 E(p) 2 r(T0(p)) 2 | T=0 


V 4 p 4 sweep 

£ 4 / 3 p 8 / 3 strain 


( 12 ) 


(13) 


where 

11 = ^ r(C) l?=0 = 2r(4) (°) r (°) + 6[r (2) (0)] 2 (14) 

is a universal inertial range constant. 

The integrals in Eq. (13) will be evaluated assuming a truncated Kolmogorov inertial 
range 

E(p) = { C k£ 2/3 P~ 5/3 for p > k 0 
\ 0 for p < ko 

Both integrals exhibit power law divergence in the limit ko —■ > 0: the integral for the 
sweeping hypothesis diverges as k 0 ; and the integral for the straining hypothesis diverges 
as k 0 . This divergence indicates as expected that sound generation is dominated by 
contributions from the most energetic scales. Defining the inertial range energy 

POO q 

K = / E(k) dk = -C K £ 2/3 k~ 2/3 

■h 0 2 

where Ck is the Kolmogorov constant, we note that in the infinite Reynolds number limit 
in which ko —■ ► 0, 

p ^ / Kl/2 sweep 
l K 5 / 2 strain 

Using these results, evaluation of the integrals in Eq. (9) leads to 


/«\5 f a sweep (V / u) 4 sweep 
\ ^strain strain 


(15) 


For the proportionality constants, we obtain 


a , 


sweep 


^strain — 


\ncV 2 

5 


'K 


9 


- 1/2 


(16) 
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The inertial range constant 1Z may be different for the straining and sweeping hypotheses, 
but this fact has not been made explicit to simplify the notation. 

The acoustic power given by Eq. (15) for strain evidently coincides with Proudman’s 
result Eq. (1). Lilley’s formula Eq. (2) also reduces to Proudman’s if the characteristic 
frequency 0 is an inertial range property, for example if 0 ~ e/K. But if 0 is chosen to 
be a sweeping frequency, say 0 ~ V/L } then Lilley’s formula coincides instead with Eq. 
(15) for sweep. 

III. Evaluation of the Proudman constant 

To evaluate the proportionality constant entrain, we use Kaneda’s 13 result 

^(0 = ex p(-°- 8l7r C 2 / 4 ) (17) 

derived for a Lagrangian modification of the direct interaction approximation in which 
time correlations are dominated by local strain. In this case, 1Z = 12.7 and therefore 

^strain = 3.6 


To evaluate the proportionality constant a swe ep, we use Kraichnan’s result 14 that the 
time correlation function for the sweeping hypothesis is 

■hm 


K0 = 

therefore 1Z = 10. Substituting in Eq. (16), 


(18) 


a !TOp = 6 C,][ 2 = 12.0 


(19) 


Woodruff has noted 15 that since the time correlation functions of Eqs. (17) and (18) 
only enter Eq. (13) through their derivatives at zero, formulas for the constant 1Z could 
be obtained from Taylor series solutions of the Navier Stokes equations. By expressing 
these derivatives in terms of single time moments without invoking quasi-normality, a 
generalization of the Lighthill 4 relation could be derived. 
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To assess the agreement between the present theory and Lilley’s, it is necessary to 
compare values of a' (QL / u) 4 and a S weep(V / u) 4 . Table I contains representative values of 
the first quantity; the Strouhal number is the factor St = QL/u of Lilley’s theory. St is 
probably 7 somewhat larger than 1.0. The value St = 1.25 which appears in Table I is an 
upper bound, not an experimental or theoretical value. Refs. 3 and 5 attempt to extract 
the parameters of Lilley’s theory from isotropic turbulence data. In Ref. 6, the sound 
radiated by isotropic turbulence was determined by numerical simulation; however, recent 
reconsideration 16 of this calculation suggests that the result may not be entirely reliable. 
With St in the range 1.00-1.25, a' St 4 takes values between 0.7 and 8.8. 

Like the Strouhal number, the ratio V/u is essentially nonuniversal since it depends 
on properties of the energy containing range. Table II contains data taken from a recent 
survey 17 of experiments with well defined inertial and energy containing ranges. The 
experiments can be identified from Ref. 17 through the value of Re\. The sweeping 
velocity V is evaluated from the energy in the scales below the 5/3 range. The values 
of 12 (V/u) 4 in Table II are generally within the range found in Table I. However, the 
data of Table I matches the entire energy spectrum including the energy containing range, 
whereas the present theory ignores the contribution of these scales to sound radiation. This 
suggests that the theoretically computed value of a sweep in Eq. (19) might be somewhat 
too large. 7 In applying experimental results 5,17 to an isotropic theory, we assume that the 
overall energetics of possibly inhomogeneous and anisotropic turbulence does not differ 
greatly from ideal isotropic turbulence. 

In both the present theory and Lilley’s, the fourth power dependence on the external 
parameters O and V is striking. In particular, the present theory predicts that decorre- 
lation by the energy containing range strongly enhances the radiation of sound; as noted 
previously, the sound radiated by the energy containing scales themselves has not been 
included in the present calculation. 
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IV. The acoustic power spectrum 


Taking the Fourier transform of Eq. (12) leads to the formula for the acoustic power 
spectrum 

poo 

p{to) ~ — / dp p~ 2 E(p) 2 r(() * r(C) lc = w /0 # _1 (20) 

c Jo 

where r(() = r(u / 6(p)) is the Fourier transform of the time correlation function defined 
by Eq. (10) and * denotes convolution with respect to the variable (. By introducing the 
change of variables q = u / 6(p) in Eq. (20), we find 

X P£ f ei/syis/s^-r/s sweep 
c 5 ( £ 5 / 2 cu -7 ! 2 strain 

The considerable difference in the scaling exponents might make an experimental test of 
the sweeping and straining hypotheses possible. The spectral scaling for sweep in Eq. (21) 
can be compared to the u~ 2 scaling 3 which has been found to give a good description 
of jet mixing noise. 18 It should be noted that Proudman’s theory 1 does not provide an 
expression for the acoustic power spectrum. 

V. Conclusions 

The calculation of sound radiation by isotropic turbulence has been considered from 
the viewpoint of analytical theories of turbulence. The calculation is accordingly based 
on universal inertial range descriptors of turbulence; empirical space and time correlation 
functions are not used. The Proudman constant is expressed theoretically in terms of the 
Kolmogorov constant, a time correlation constant, and various purely geometric factors. 
The calculation shows a significant effect of the behavior of Eulerian time correlations 
on sound radiation. If decay of time correlations is dominated by local straining, then 
sound radiation is given by Proudman’s original formula 1 and the acoustic power spectrum 
follows an cu -7 / 2 law; whereas the more plausible hypothesis that decay of time correlations 
is dominated by sweeping by the largest energetic scales results in a formula much closer 
to Lilley’s 3 , and the acoustic power spectrum follows an cu -4 / 3 law. 
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TABLE I 


a' St a' St 41 
Ref. 3 3.2 - 3.6 1.25 7.8 - 8.8 
Ref. 5 0.7 - 2.7 1.25 1.7 - 6.5 
Ref. 6 2.6 


TABLE II 


Re x 

V/u 

12 (V/u) 

3180 

.48 

0.6 

1450 

.86 

6.4 

850 

.65 

2.2 

540 

.59 

1.5 

308 

.84 

5.9 
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